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Abstract 

Supersymmetric mass spectra within two variants of the seesaw mechanism, commonly known 
as type-II and type-III seesaw, are calculated using full 2-loop RGEs and minimal Supergravity 
boundary conditions. The type-II seesaw is realized using one pair of 15 and 15 superfields, while 
the type-III is realized using three copies of 24^ superfields. Using published, estimated errors on 
SUSY mass observables attainable at the LHC and in a combined LHC-I-ILC analysis, we calculate 
expected errors for the parameters of the models, most notably the seesaw scale. If SUSY particles 
are within the reach of the ILC, pure mSugra can be distinguished from mSugra plus type-II or 
type-III seesaw for nearly all relevant values of the seesaw scale. Even in the case when only the 
much less accurate LHC measurements are used, we find that indications for the seesaw can be 
found in favourable parts of the parameter space. Since our conclusions crucially depend on the 
reliability of the theoretically forecasted error bars, we discuss in some detail the accuracies which 
need to be achieved for the most important LHC and ILC observables before an analysis, such as 
the one presented here, can find any hints for type-II or type-III seesaw in SUSY spectra. 
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I. INTRODUCTION 



In the minimal supersymmetric extension of the standard model (MSSM) all soft SUSY 
breaking mass terms are treated as free parameters, to be fixed at the electro-weak scale. 
However, these soft parameters potentially contain a wealth of information about physics 
at the high scale and understanding the nature of SUSY breaking will become the main 
challenge, if signals of SUSY are found at the LHC Highly precise mass measurements 
will be needed to distinguish between different SUSY breaking schemes such as "minimal 
supergravity" (mSugra) [H, |2|, anomaly mediated SUSY breaking (AMSB) js], 0] or gauge 
mediated SUSY breaking (GMSB) [sl, to name just the most familiar ones. 

However, all of the models mentioned above break SUSY at energies inaccessible for 
collider experiments. Thus, theoretical extrapolations from the TeV scale to the high energy 
scale will be needed and any "test" of SUSY breaking schemes can at best take the form 



of a consistency check. Based on the results of [6H8| detailed calculations have been done, 
quantifying the accuracy with which such tests can be done using data from LHC and a 
possible ILC j9-13|. However, these works concentrated on models with MSSM particle 
content and thus did not attempt to take into account the observed non-zero neutrino 
masses. In this paper we study the prospects for the LHC and for a combined LHC+ILC 
analysis for finding indirect hints for the presence of a high-scale seesaw mechanism in SUSY 
spectra. 

The MSSM assumes that R-parity is conserved and thus, jus t as in the standard model, 
neutrino masses vanish. Neutrino oscillation experiments 
that at least two neutrino masses are non-zero 
neutrino masses, "the seesaw" mechanism 19 



14 



17 



however, have shown 
Among the myriad of possible models of 
21| is undoubtedly the most popular one. The 
classical version of the seesaw 19|, |20j introduces (at least two) fermionic singlets ( "right- 
handed neutrinos") with some large, but arbitrary Majorana mass Mr. The smallness of 
the observed neutrino masses is then a straightforward consequence of being inversely 
proportional to Mr. This variant of the seesaw is now usually called type-I seesaw. 

At tree- level there are only three realizations of the seesaw mechanism 22|]. In addition 



24 



to type-I, the seesaw can be generated by the exchange of a scalar triplet (type-II) [23 
or by a fermionic triplet in the adjoint representation (type-Ill) 25|. Common to all of 
them is that for ~ -^/AmJ ~ 0.05 eV, where Am^ is the atmospheric neutrino mass 
splitting, and couplings of order 0(1) the scale of the seesaw is estimated to be very roughly 
mss ~ 10^^ GeV. 

Extending the standard model (SM) with a seesaw mechanism leaves no experimental 
signal apart from the neutrino masses themselves. The situation is different in the super- 
symmetric seesaw. There are two kind of measurements which potentially can give indirect 
information about the seesaw parameters: Lepton flavour violating (LFV) decays and su- 
perpartner mass measurements. 



The literature on LFV in SUSY seesaw is vast 26[. It was pointed out already in 27 



that LFV is practically unavoidable in SUSY seesaw, even if the SUSY breaking boundary 
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conditions are completely flavour blind. However, for (3 generation) type-I and type-Ill 
seesaws the seesaw mechanism has more free parameters than there are observables in the 
neutral and charged lepton sectors. The three LFV entries in the left slepton mass matrix 
can then be made arbitrarily small (or large) independent from any neutrino physics and 
thus there are no definite predictions for LFV decays in SUSY seesaw. ^ Indeed, if any 
charged LFV is ever observed one would probably turn the argument around and try to 



33-35 



learn indirectly about the unknown seesaw parameters instead 

In contrast, there are only very few papers, which have studied the impact of the seesaw 
on SUSY particle masses. Some aspects of type-I seesaw have been studied focusing on what 



can be learned from precision measurements in the slepton sector 10|, I36l-l38|. Moreover, in 
such a scenario a splitting between the masses of the selectrons and smuons can occur which 



might be measurable at the LHC 



39 
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Changes in SUSY spectra can lead to changes in 



the expected relic density for the cold dark matter. The impact of large values of soft terms 
in the sneutrino sector 
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and of large values for the trilinear Aq parameter [43|, |4J] have 
been studied in this context. 

The relative scarcity of publications on SUSY spectra and the seesaw is probably ex- 
plained by the fact that type-I seesaw, the undoubtedly most popular variant, adds only 
singlets to the MSSM particle content. If the Yukawa couplings of these singlets are smaller 
than, say, the gauge couplings any effects of the right-handed neutrinos on the SUSY mass 
eigenvalues become negligibly small. This leaves only a rather small window for the seesaw 
scale, mss, say, roughly [4 x 10^^, 1.2 x 10^^] GeV where any measurable shifts in SUSY 
masses can be expected at all. And it is, of course, exactly this range for 77155 where the 
largest values for LFV decays are expected. Exceptions from this general rule can be found 
in models where one departs from the universality assumption of the mSUGRA parameters 



with huge soft SUSY breaking parameters in the seesaw sector one gets larger effects [45 
in particular in the Higgs sector 
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Changes in SUSY spectra with respect to, say, mSugra expectations are expected to be 
much larger in type-II and type-Ill seesaws, but very little work has been done on these 
seesaw variants as well. Type-II and DM has been studied in 47l|, while for a study of 



the type-Ill seesaw with emphasis on spectra and LFV see [481]. In 49|] it was pointed 
out, that one can form different combinations of soft SUSY breaking parameters, which 
at 1-loop order do not depend on the mSugra parameters. Consistent departures of these 



invariants" from mSugra expectations could then be taken as indirect hints of the seesaw. 



(See, however, [29l. l48l| for the importance of 2-loop effects on the "invariants".) The above 



papers [29|, |48|, |49| have pointed out, how type-II and type-Ill leave traces in SUSY spectra, 
in principle. They did not, however, attempt to quantify the accuracies needed to find hints 
of the seesaw in experimental mass measurements. To our knowledge the current paper is 



^ The situation is different in "minimal" typc-II seesaw. Here, ratios of different LFV decays are related 
to neutrino angles, if (a) mSugra boundary conditions are assumed 
seesaw triplet is also responsible for SUSY breaking, see 30l432|. 
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2^ or (b) in schemes where the 
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the first step in this direction. 

Here we calculate the low-energy SUSY spectra for type-II and type-Ill seesaw in mSugra 
and confront our theoretical results with expectations for the accuracy of SUSY mass mea- 
surements at the LHC and at a possible combined LHC-I-ILC analysis 0, S]- Given the 
estimated errors on SUSY masses obtained in detailed simulations 0, g we calculate ex- 
pected x^-distributions for the two different seesaw models, in order to give a theoretical 
forecast on the expected errors on the model parameters, most notably the error on the 
"determination" of the seesaw scale mss- 

Our main results are the following: With the highly accurate mass measurements ex- 
pected for the ILC it should be possible to distinguish between pure mSugra (i.e. a model 
with no seesaw at all) and mSugra plus type-II or type-Ill seesaw for almost any relevant 
value of mss, if at least some SUSY particles are within kinematical reach of the ILC. We 
find it noteworthy, however, that even with the much less accurate data, expected from 
LHC measurements only, it seems possible to distinguish pure mSugra from the mSugra 
plus seesaw models in some favorable parts of parameter space. Obviously, all our results 
depend crucially on the - currently only theoretically estimated - errors, with which SUSY 
mass can be measured at LHC and ILC. We therefore discuss in some details what are the 
observables needed and the required error bars on these observables, before an analysis, such 
as the one presented here, can find any hints of the seesaw in SUSY spectra. 

The rest of this paper is organized as follows. In the next section we summarize the three 
different variants of the seesaw, to set up the notation. We embed the new particles required 
by the different seesaw mechanisms in complete SU{5) multiplets in order to maintain the 
successful unification of gauge couplings observed in the MSSM. Section 3 contains the bulk 
of this paper. It first defines our setup, lists the observables we use and then discusses SUSY 
spectra in the different models. With these results we then proceed to calculate theoretical 

distributions. We first discuss a combined LHC+ILC analysis and then go to the case 
of using only the less accurate LHC data. Our results, of course, depend crucially on the 
accuracy with which the SUSY masses can be measured in future accelerator experiments. 
We therefore dedicate a subsection to discuss in detail the requirements for the accuracies 
on the most important observables needed for our analysis. We then close with a short 
discussion and outlook. 



II. SUPERSYMMETRIC SEESAWS 

In this section we briefly recall the main features of the three tree-level variants of the 



seesaw. A more detailed discussion including the embedding in SU{5) can be found in [50 
For brevity, we will discuss only the superpotential terms. In all cases, we start with the 
MSSM superpotential 



Wmssm = U'YuQ ■ - D'^YaQ ■ Ha - E'^Y.L ■ Ha + fiH,, ■ Ha . (1) 
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Here, A- B = A1B2 — A2B1 denotes the SU{2) invariant product of two SU{2) doublets and 
Ya are (3, 3) Yukawa coupling matrices, as usual. Only for completeness we mention that a 
seesaw type-I is obtained by adding to the above superpotential: 

Wj = N'Y.L ■ if„ + ]^N'MrN'' . (2) 

Integrating out the heavy singlets, after electro-weak symmetry breaking, eq. ([2]) leads to 
the famous seesaw formula 



V. 
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2 YjM-'Y,. (3) 

The atmospheric neutrino mass scale requires at least one neutrino to be heavier than 
Am\ ~ 0.05 eV. For Yukawas of order 0{1) this results in a "seesaw scale" of roughly 
Mr ~ 10^^ GeV. Since eq. ([2]) only adds singlets to eq. ([1]) one does in general not 
expect any sizeable changes in the supersymmetric mass spectra at the electro- weak scale, 
see also the discussion in the introduction and in 48|. In section UlIB I we will discuss, how 
the introduction of non-singlets changes the expected SUSY spectra and comment on the 
expected change for type-I. 

In supersymmetric models the simplest way to generate a type-II, while maintaining 
gauge coupling unification, is to add a pair of 15-plets of SU{5) to eq. ([1]). The SU{5) 
invariant superpotential than reads 

W = ^Yi55m ■ 15 ■ 5a/ + ■15-5h + ^^25^ ■ 15 ■ 5^ + Y5IO ■ 5 ■ 5^^ 

+ YioIOm ■ 10m ■ 5h + M15I5 ■ 15 + M,5h ■ 5^- (4) 

Here, 5*/ and 10^/ are the usual SU{5) matter multiplets and 5h = {H'^,Hu) and 5h 
{H\Hd). Under SU{'i) x SU{2) x U{1) the 15-plet decomposes as 



15 = 5 + T + Z , (5) 

S ~ (6,1,-^), r~(l,3,l), Z~(3,2,^). 

3 D 

Below the GUT scale, M^, in the 5't/(5)-broken phase the superpotential reads 

Wii = -^{YtLTiL + YsD''SiD^) + YzD%L 
V2 

+ -^(Aii/,fiif, + \2H^%H^) + M^fifs + MzZ^Z2 + MgSA. (6) 

where fields with index 1 (2) originate from the 15-plet (15-plet). The first term in eq. (E]) 
is responsible for the generation of neutrino masses, which at low energies are given by 

Similar to type-I, the seesaw scale is estimated to be ^ ^ lO^^GeV. 
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In the case of a seesaw model type-III one needs new fermions E at the high scale belonging 
to the adjoint representation of SU{2). The simplest complete SU{5) embedding possible 
is the 24-plet 22|. The superpotential of the unbroken SU{5) is then 

W = V25mY^10m5h - ^10mYio10m5h + 5h2AmY^''5m + UamM2^2Am • (8) 
The 24m decomposes under 5f/(3) x SU{2) x f/(l) as 

24m = (1, 1, 0) + (8, 1, 0) + (1, 3, 0) + (3, 2, -5/6) + (3*, 2, 5/6) , (9) 
= Bm + Gm + Wm + Xm + Xm ■ 

The Bm has the same quantum numbers as the iV^, while the fermionic component of the 
Wm corresponds to the S. Thus, the 24m always produces a combination of the type-I and 
type-III seesaw. 

In the SU (5) broken phase the superpotential contains 

Win = Hu{WmYw - ^BmYb)L + H^XmYxD' 

+ ^BmMbBm + ^GmMgGm + ^WmMwWm + XmMxXm- (10) 
Integrating out the heavy fields, as before, leads to 



m 



-| {to^^^^'^' + lywM^'yw) ■ (11) 

There are two contributions: (i) from the gauge singlet and (ii) from the SU(2) triplet. 
Starting with a common Y^^^ , Y^ evolves slightly differently from Y^ under the RGEs. 
Thus, in principle two non-zero neutrino masses are generated from one 24m only. However, 
the ratio of the two non-zero neutrino masses generated in the RGE running is much too 
tiny to explain the observed neutrino data and thus at least 2 copies of 24m are needed for 
a realistic neutrino mass spectrum. In our numerical calculations we use 3 copies of 24m, 
motivated by the observed 3 generations. With VY^ /Y^ ~ 1, one can simplify eq. ( ITTl) to 

mu = -vl^Y^M^'Yw (12) 
The scale of Mw is then estimated to be mss ~ 8 x 10^^ GeV for = 0(1). 



III. RESULTS AND DISCUSSION 

In subsection IIII Al we will define our setup and discuss the input observables. In IIII Bl 
we discuss the SUSY spectra in the different models and how the observables change under 
changes of the seesaw scale. Section [III CI we present our results for a combined LHC+ILC 
analysis, while IIII Dl shows the results for an analysis using only LHC data. Section IIII El 
discusses the accuracies on the different observables, which need to be achieved experimen- 
tally, before any conclusions on the presence (or absence) of a seesaw mechanism can be 
drawn. Given the inherent unreliability of theoretical error forecasts in general, section IIII El 
can be considered to contain the central parts of the current paper. 
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A. Setup, observables and data input 



All the numerical results shown in the following have been obtained with the programme 



package SPheno 5l|, |52|. The RGE equations, complete at the 2-loop order, have been 



calculated and incorporated into SPheno with the help of SARAH 53|-|55|. Details and 



discussion of the implementation can be found in |48 . 

To completely specify the low-energy SUSY spectra, we have to assume a specific SUSY 
breaking scenario. In this paper we use mSugra. mSugra is defined at the GUT-scale, Mq, 
by: a common gaugino mass M1/2, a common scalar mass mo and the trilinear coupling 
Ao, which gets multiplied by the corresponding Yukawa couplings to obtain the trilinear 
couplings in the soft SUSY breaking Lagrangian. In addition, at the electro-weak scale, 
tan/3 = Vu/vd is fixed. Here, as usual, and u„ are the vacuum expectation values (vevs) 
of the neutral component of and Hu-, respectively. Finally, the sign of the /i parameter 
has to be chosen. 

In the following we will call "pure mSugra", pmSugra for short, the version of the model 
with no seesaw mechanism at all. Note that this is equivalent to putting the seesaw scale 
mss equal to the GUT scale Mq- For brevity, we will call "tpye-II" and "type-Ill" mod- 
els with mSugra boundary conditions, to which on top of the MSSM particle content the 
corresponding "seesaw particles", as specified in the previous section, are added at scale 
mss- 

In our numerical calculations we concentrate on some selected sets of mSugra parameters. 
This is mainly motivated by the exorbitant amount of CPU time a full scan over the mSugra 
space would require, see also the discussion below. The points we have studied are SPSla' 



[81] and the points SPSlb and SPSS [56|]. In addition, for reasons explained in section 
IIIIBl we consider a few more points with modified mSugra parameters. We will call these 
points MSP-1 (mo,Mi/2,tan/3,Ao) = (70, 400, 10, -300), MSP-2 (220,700,30,0) and MSP-3 
(120,720, 10,0). MSP-1 is similar to SPSla' but with a larger value of M1/2, MSP-2 is a 
point with larger tan /3 and MSP-3 is similar to SPS3, but again with a larger value of Mi/2- 
All these points choose /i > 0. We have not found any qualitatively new features in points 
with negative /i, as far as the determination of mss is concerned. 

Observables and their theoretically forecasted errors are taken from the tables (5.13) 
and (5.14) of [3] and from j^. For the LHC we take into account the "edge variables": 
[muf^^e^ {miq)l^ {miqf^^l, ("^iig)edge aud (m;ig)thresh from the decay chain ql -> xIq and 



X2 ^ II ^ llXi |57H59| . In addition, we consider {min,)thresh, ("^r+r-) (from decays involving 
the lighter stau) and the mass differences A-^, = mg — mf^,, with i = 1,2, A^^^o = ~m^Q 
and A;-^^o = m^-^ —my.o. Since m^^ ^ m^^ and rriun — m^^ ^ m^^ ~ m^^ applies for a large 
range of the parameter space LHC measurements will not be able to distinguish between the 
first two generation sfermions. ^ This allows us to define the masses mj^ = {m^^ + m^^)/2 
and m^^ = (m^^^ -|- mj^ + m^^ + ttLsj^)/^ which will be used from now on for the mass 



^ See however |39l441|. 
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differences A^^^o and A^-^^o. As discussed below in section UlI Et especially A-^^ and the 
edge variables are important for the LHC analysis. For the ILC we assume that at least 
m^o, rrien — and mg^ ~ m^^ are kinematically accessible. In addition, whenever within 
the reach of the ILC, we also take into account fi, xt ^^nd ii, which are, however, less 
important. We also assume that the lighter Higgs, has been found and its mass measured 
with an accuracy which depends on whether the analysis is for LHC only or for LHC+ILC, 
see the corresponding error estimates in Errors for the ILC are taken directly from 
the tables of the above papers. For the error bars for the LHC, however, we have rescaled 
all statistical errors from the values for 300 fb^^ to a luminosity of (only) 100 fb^^. To 
be conservative the total error is obtained summing statistic and systematic error linearly. 
Note that we did not make use of the combined LHC and ILC errors calculated in the papers 
mentioned above. When we discuss the calculations for LHC and ILC observables in llll CI we 
refer to an analysis in which the LHC and ILC observables are all enabled but the errors are 
the errors for the LHC or ILC only. Nevertheless, we have checked that using the combined 
errors changes the results only by an irrelevant amount. We will call this analysis therefore 
"ILC+LHC combined". 

In [3| and [sl, only standard SPS points have been studied in detail. In the calculation 
of the x^-distributions we assume that relative errors for different mSugra points and/or 
seesaw points are constant. The assumption to use constant relative errors in all of our 
calculations is, of course, a crucial simplification which has to be checked very carefully. 
However, we have chosen to do so for the following two reasons: (a) It allows us to perform 
a analysis for all different spectra within reasonable CPU time. And (b) uncertainties of 
the theoretically forecasted error bars are nearly impossible to estimate. Only experiments 
can finally determine total errors on observables. We thus use errors-as-predicted and discuss 
in section IIIIEt how our conclusions will change as a function of these unknown errors. 

To numerically estimate the allowed ranges for the model parameters we use a simple 
procedure. We have found that, see below, errors on mo, M1/2 and niss are very strongly 
correlated. To assure that our estimates are reliable in all cases we have written two com- 
pletely independent numerical codes. The first of these is based on MINUIT, ^. enforcing 
the russ scan while MINUIT is fitting the parameters mo, M1/2, tan/3 and Aq for fixed mss- 
The second code uses a straight-forward but slow Monte Carlo random walk procedure, 
which can be "heated" to find separated minima. In the MC calculations we use usually 
a (few) 10^ points to assure convergence. This makes the MC code slow, but reliable. We 
have done calculations using both codes in all cases necessary, to ensure that convergence 
has been reached. 

Finally we need to mention that in all calculations shown below we put neutrino Yukawa 
couplings to negligibly small values, unless noted otherwise. Again the reason for this choice 



^ Minimization package from the CERN Program Library. Doeumentation can be found at 
http : / / cernlib . web . cern . ch/cernlib/ 
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is simply to limit the amount of CPU time necessary for our fits. ^ We will comment, 
however, in section flIIIE|) on the differences expected with fits, where the Yukawas are 
chosen to fit neutrino data correctly. In general, iimss is below, say, 10^^ GeV the differences 
of the full fit to our calculation with negligible Yukawas is found to be completely irrelevant. 
If any hints of seesaw with mss in the range [10^^, 10^^] were indeed found in SUSY spectra, 
however, we expect that a full analysis would find results which differ by some (10 — 30) % 
(depending on the exact value of mss) from our preliminary numbers. 



B. Mass spectra and LHC and ILC observables 



In this subsection we briefly summarize the differences in the calculated mass spectra of 
the different models and how this affects the LHC and ILC observables, which we use in 
our fits. In the following all numerical results shown in the figures have been calculated 
solving the full 2-loop RGEs numerically. Moreover, we have taken into account the 1- 
loop thresholds of the seesaw particles at the seesaw scale as described in |48| and include 
the one-loop contributions to the SUSY masses 
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We have also included the shifts of 
the gauge couplings to the DR-scheme. It is instructive to discuss some (semi-) analytical 
approximations at 1-loop order, which will allow to understand qualitatively the numerical 
results. We stress that none of the following approximations is used in any way in the 
numerical calculations. 

The introduction of complete SU{5) multiplets at a scale below the GUT scale changes 
the running of the gauge couplings. At 1-loop order the gauge couplings at the different 



scales are given as 29 
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-7) for SM and b = (6i, 62,^3 
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for MSSM. mss denotes the seesaw scale, i.e. the mass of the 15-plet or the mass(es)^ of 



^ We let 5 parameters flow freely. For a type-II, for example, we have in Yt six more complex parameters. 

A full fit would require minimizing for 5+12-3=14 parameters, which can not done for all spectra we 

need to consider within realistic amounts of CPU time. 
^ We assume that the 3 copies of 24-plets are degenerate. In principle, given enough accurately measured 

observables, it might be possible to drop this assumption. 
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the 24-plets. For the case of the 15-plet one finds Abi = 7 whereas for the case with three 
24-plets one finds Abi = 15, since each 24m gives a Abi = 5. It is easy to show that at 
the 1-loop level the GUT scale is not changed by the introduction of the complete SU{5) 
multiplets. However, the Abi 7^ lead to a faster "running" of the gauge couplings and thus 
to a larger value of ^(Mg) compared to the mSugra case. For seesaw scales smaller than 
roughly niss ~ 10^ GeV (10^^ GeV) one then encounters a Landau pol in a^Mc) for seesaw 



type-II (type-Ill) 
on perturbativity. 
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This defines in each case a lower limit on the seesaw scale, if we insist 



Gaugino masses evolve like gauge couplings: 



MiimsusY) 



a{Mc 



Ml/2. 



(14) 



Eq. f|T^ implies that the ratio M2/M1, which is measured at low-energies, has the usual 
mSugra value, but the relationship to M1/2 is changed. I.e. since a^Mc) is larger in the 
seesaw case than in the standard pmSugra, Mj are smaller in seesaw than in pmSugra. 
For the soft mass parameters of the first two generations one obtains j29j 
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The various coefficients cf are given in table [B 




/2' 



(15) 



Abi) log 



Ml 



m 



ss. 



(16) 



/ 


E 


L 


D U Q 




6 


3 


2 8 1 


^\ 


5 


10 


15 15 30 


J 





3 
2 


1 


'^3 








coloo 
coloo 
ooloo 



TABLE I: Coefficients c{ for eq. ([I5 



In the limit mss — ^ Mq the functions fl go to zero and one recovers the standard mSugra 
estimations for the sfermion masses. For any rriss below Mq, the contribution from /j are 
smaller than in the mSugra case, due to the prefactor which is always smaller than one. The 
contribution from the // can only partially compensate for this and thus, at low energies 
for a given pair of mo and M1/2 one expects the sfermion masses to be smaller in seesaw 



^ Note, however, the SPheno never allows us to push the seesaw scales down to these limits. Convergence 
problems are usually encountered already for q;(A/g) > 0.25 depending on the mSugra parameters. 
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FIG. 1: Supersymmetric masses for two specific choices of mSugra parameters as a function of the 
seesaw scale. To the left: type-II. To the right: type-Ill. The mSugra parameters are fixed at the 
values of SPSla' (top) and to MSP-1 (bottom). Note the different scales for type-II and type-Ill. 



tlian in pmSugra. It is important, however, to note that the different coefficients q differ 
not only from sparticle to sparticle, but also are different for the same particle but different 
gauge groups. This observation is fundamentally the reason explaining the statement that 
accurate sfermion mass measurements will allow to distinguish pure mSugra from mSugra 
plus seesaw.^ As in the case of pmSugra, coloured particles are expected to be heavier than 
non-coloured ones in seesaw. 

Before discussing the numerical results, we briefly comment on seesaw type-I. In type- 
I one only adds singlets to the MSSM particle content. Thus, Abi = and there is 
no deformation of the spectrum with respect to mSugra due to the gauge part. The only 
change one expects for type-I is due to a different running of m|, when Yukawas are taken 
into account. Since only m| is affected, most of the observables we have discussed above 
are not sensitive to the seesaw scale in type-I and our current analysis can not directly be 
applied to type-I seesaw. 

Fig. ([1]) shows some examples of SUSY masses for two specific choices of mSugra pa- 



If all Ci where the same, one could always fit the data by a simple rescaling of toq ^-nd Mi/2- 
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FIG. 2: Running masses as a function of the seesaw scale, left: type-II; right: type-Ill. As fig. 
p!|). but showing only the point SPSla'. The masses shown in this figure are the most important 
coloured particles for our analysis. 



rameters. The top panel shows mSugra parameters chosen as in the standard point SPSla' 
[sl, while the bottom panel shows MSP-1 for comparison. All plots show masses as function 
of mss, to the left for seesaw type-II and to the right for type-Ill. Shown are masses of 
the lighter two neutralinos, Xi and X25 masses of charged sleptons. As is also typical 
in mSugra, the mass of the lighter chargino is very similar to m^o and smuons and 
selectrons are nearly degenerate. We note in passing, that the mass of the lightest Higgs, 
h'^, shows little or no sensitivity at all on 77155 but is important in the fits. 

As discussed above, all masses get smaller for smaller values of 77755 and are always smaller 
than in the pmSugra limit. Note the wide range for type-II shown and the much smaller 
range of 77755 plotted for type-Ill. Ratios of gaugino masses follow standard expectations for 
all values of 77755 and both types of seesaw. The slopes of the curves are different for different 
sparticles and the relative changes are larger in type-Ill than in type-II. This simply reflects 
the fact that type-Ill causes a larger change in the beta coefficients (A6j = 15) than type-II 
{Ah = 7). 

For the point SPSla' the lighter chargino becomes lighter than 105 GeV for type-II (type- 
Ill) seesaw scales below roughly 77755 ~ 2 x 10^^ GeV {mss ~ 10^^ GeV). Thus, 77755 below 



these values are ruled out by the LEP bounds 6l|, l62| . Note that this implies that type-Ill 



can not explain neutrino data for mSugra parameters as in SPSla' with Yukawas smaller 
than < 1. 

Changing the seesaw scale can lead to a different mass ordering for different sparticles. 
For example, for SPSla' the X2 is heavier than and fi in the mSugra limit (seesaw scale 
equal to Mq), but lighter than for type-II (type-Ill) seesaw scales below 77755 ~ 3 x 10^^ 
(77755 ~ 8 X 10^^). This is important for our study, since as a function of 77755 it can happen 
that some observables are kinematically open for some values of 77755 but not for others, see 
also the discussion in section IIIIDi 

The modified value of M1/2 = 400 GeV in MSP-1 with respect to SPSla' is motivated 
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by the fact that for this choice of parameters the edge variables from the chain X2 ~^ ^rJ- ^ 
llXi are kinematically possible for all relevant values of rriss- The larger value of M1/2 
implies heavier neutralinos and also that the LEP bounds on sparticle masses are fulfilled 
for all values of ^5-5 shown. Note, that all sparticles shown in the plot are kinematically 
accessible at an ILC with ^/s = 1 TeV. For MSP-1 the lighter stau is the LSP for mss larger 
than roughly mss = 10^^ GeV. Thus, this point formally has no cosmologically acceptable 
pmSugra limit. 

Fig. ([2]) shows the dependence of several coloured sparticle masses on mss- Again to 
the left (right) we show seesaw type-II (type-Ill). Note the different scales for type-II and 
type-Ill. We show only the values for SPSla' in this figure, masses for MSP-1 are larger 
but behave qualitatively very similar. The relative change of masses as a function of mss 
is much larger than for the non-coloured sparticle masses shown in fig. ([1]). Here the range 



where ~ m^^ ^ < 300 GeV is excluded by Tevatron data [6II, 163|. However, this region 
is already excluded by LEP data. Note that mg > m^^ for all values of mss in this point. 
Coloured sparticle production gives the bulk of the SUSY cross section at the LHC as usual. 
In these points most of the coloured sparticles are not kinematically accessible at the ILC, 
except for low values of mss- Except ti we therefore do not take into account measurements 
of coloured sparticles at the ILC in our analysis, even though they could be potentially much 
more accurate than the corresponding measurements at the LHC. 

With the masses shown in fig. ([T]) and fig. ([2]) one obtains the LHC observables shown in 
fig- (E]) for type-II. Again in the top panel we show SPSla' and in the bottom panel MSP-1. 
The figure shows several mass differences (left) and the edge variables (right) stemming from 



the decay chain q — )■ with the subsequent decay X2 ~^ ^^^^ ~^ ^^^^X? |57H59|. We have 
normalized all observables to their expected values for mss = Mq- Thus relative changes in 
the different observables with respect to pmSugra are plotted. 



The two kinks in the running of (mig)'^^^'^ and (mig)'^^^ stem from the fact that one has 
to consider different cases in these observables. They can be written as 58 

= max[(m— ^)^ (m^ )2] (17) 

2 2 ^ 
ip: — m 

1 X2' 



{mi,)t^ = min[(m;:-J^K-<o)(m? -mJo)/(2m? -mjo)] 



where 



imZZ,r = H - ^Jo)(mJo - m?J/mJo (18) 

= K--Jo)K--Jo)/(mfj. 

These conditions change as a function of mss causing the kinks shown in the figure. Except 
for the {muY^^^ different cases appear in the expressions for all edges, but only for the 
variables {miq)l^^ and (^«ig)hfgh do the kinematical conditions change as function of 777,55 
normally. 

The plot in fig. ([3]) demonstrates the strong dependence of the LHC observables on mss- 
Increasing and decreasing values of the edges are possible, while mass differences usually 
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FIG. 3: Relative change of different LHC observables as a function of the seesaw scale for type- 
II seesaw. Top mSugra parameters as in SPSla', bottom MSP-1. For an explanation see text. 
(mee)'^'^^° is repeated in the left plot for comparison. 



decrease for lower values of 77255. Note that in the fits, discussed in the next subsections, 
observables which show (a) the largest relative change with respect to 77155 and (b) have 
the smallest expected errors will give the most important contributions. Finally we mention 
that fig. (|3]) shows only type-II, since results for type-Ill are qualitatively similar but with 
larger relative changes. 

C. analysis for combined LHC and ILC data 

In this section we take into account all possible LHC and ILC observables. We discuss this 
more futuristic (but simpler) case first. Results for an analysis taking only LHC observables 
are discussed in the next subsection. We note in passing that we have checked that we can 
roughly reproduce the error on parameters for the pure mSugra results for the point SPSla' 
discussed in detail in jsl. 

Fig. (jl]) shows the allowed ranges of the parameters for the point MSP-1 and one specific 
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FIG. 4: Allowed parameter space for mo, M1/2, tan/3, ^0 ^^nd mss for all 5 parameters varied 
freely. The input value for the seesaw scale is mss = 5 x 10^'^ GeV and seesaw type-II has been 
used. 



choice of 777.55 = 5 x 10^^ GeV for type-II seesaw. The allowed regions have been found in a 
MC random walk procedure letting 5 parameters, mo, M1/2, ^0, tan/3 and 77755, float freely. 
The ranges shown correspond to a Ax^ — 5.89, i.e. 1 a c.l. for 5 free parameters, where 
we have taken into account the correlations between the various parameters. Plotted are 
different 2-dimensional projections of parameters. 

As mentioned already above, the three parameters ?77o, M1/2 and 77755 are highly corre- 
lated among each other. Lower values of 77755 can be compensated by increasing M1/2 and 
decreasing 7779 at the same time. This feature is present in all parameter space for both 
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FIG. 5: distributions of the random walk for MSP-1 and Seesaw type II. The dashed hne 
indicates a of 5.89. Recall that 5.89 corresponds to a la confidence level for five free parameters. 
The plots show the distributions for 777,55 = 1 x 10^^ GeV (left) and 771,55 = 1-3 x 10^^ GeV 
(right). 



types of seesaw. This correlation results in errors on mo and M1/2 which are larger (some 
times much larger, see below) than in pmSugra for the same input errors on observables. 
We note that the in this calculation is dominated by the much more accurate ILC data, 
see also the discussion in section ( IIII El) below. 

In contrast, tan/? and Aq show very little correlation with 777,55 (and ?t7o and M1/2) and 
only a rather moderate correlation among themselves, tan /3 and Aq are mostly determined 
by the Higgs mass measurement, and to some extend by 3rd generation sfermions. Note 
that Aq and tan/3 do not have much influence on the determination of rriQ, M1/2 and mss, 
apart from a slight increase in the errors of the latter. However if 7T7o cannot be fixed 
a determination of Aq and also tan /3 becomes practically impossible, because almost any 
shift of tan/3 and Aq can then be compensated by changing tuq and/or Mi/2- This will be 
important when we discuss the calculations using LHC observables only in section IIII D[ 

For this choice of parameters, the error on 77755 itself is found to be around A77755 ~ 
1.2 X 10^'^, i.e. values of 77755 = Mg are formally excluded by many standard deviations. 
However, A77755 is a very strong function of 77755 itself, as we will discuss below. 



16 



Seesaw II (mo : 70, M1/2 ■ 400, tan/3 : 10, Aq : -300) 
77.5 



Seesaw II (m,, : 70, Nh/2 : 400, tan/3 : 10, ^0 : -300) 



> 

o 




10^5 2 X 10^'"' lO^*' 2 X 10" 

mss [GeV] 





425 - 




400 - 


> 










375 - 












350 - 




325 - 




10^'^ 2 X 10^'^ lO^*^ 2 X 10^'' 

mss [GeV] 



FIG. 6: Allowed ranges of parameters rngg, uiq and Mi^2 MSP-1 and type-II seesaw with input 
77155 = 1.3 X 10^^ GeV. Two separate solutions appear, one fake but acceptable minimum is at 

mss = Mg- 



Fig. (jl]) shows results for a comparatively low value of 77755. Fig. ([5]) shows distribu- 
tions for MSP-1 obtained by a random walk for seesaw type II and two slightly different but 
much higher values of mss'- To the left: mss = 10^^ GeV and to the right: 777,55 = 1-3 x 10^^. 
The plots show the true x^-minimum and a second (fake) minimum at 77755 = Mq- For the 
lower value of 77155 this fake minimum is just excluded at 1 a c.L, while for the slightly 
higher value of 77755 it is accepted at 1 a c.l. These kind of false minima appear in all our 
calculations when 77755 approaches Mq. This is to be expected, since the models approach 
pmSugra in this limit. 

Fig. (El) shows the allowed range of parameters 77755, "^o and M1/2 for 77^55 = 1.3 x 10^^ 
GeV. Two separate minima show up. For slightly larger values of 77755 the two solutions 
overlap completely. Note that this also increases the errors on and M1/2. For slightly 
smaller values of mss this fake solution disappears resulting in a drastic decrease in the 
error bars of these three parameters. In case of type-Ill this fake minima does not show 
up separately, but indirectly by deforming the distributions. Thus for type-Ill the error 
bars go up to Mq until it gets compatible with pmSUGRA. This will be important when 
we discuss mSUGRA plus type-Ill later on. 

In case of the ILC-I-LHC analysis this kind of "false" minima are usually the only class 
of fake minima that appear. Using only LHC data, the distributions are not that well 
behaved and false minima can also appear considerably below Mg, this will be shown in the 
next subsection. 

Fig. ([7]) shows 1 cr, 3 a and 5 a c.l. error bars on the different parameters of the model 
as a function of the seesaw scale for one specific mSugra set, MSP-1, for the case of type-II. 
The plots show a large range of 77755 between [10^^, 10^^] GeV. Lower values of 77755 are 
in principle possible, but show no new features. Larger values of 77155 can not fit current 
neutrino data. Error bars on all parameters increase with increasing values of 77755, and for 
values of 77755 larger than (roughly) (1 — 2) x 10^^ GeV the error A(r?755) is so large that the 
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FIG. 7: Error of 77155, '"^O; -^1/2; to77/3 and against 77155 for all 5 parameters freely varied. For 
these plots we used the LHC and ILC observables. The chosen values for 77755, "t-Oj -^1/2? tcLufi 
and Aq are the values for MSP-1. The plots show the results for seesaw type II where we used a 
la, 3a and 5a c.l.. 



type-II can no longer be distinguished from pmSugra at the 1-a level, given the ILC+LHC 
observables with our "standard" errors. A (77755) decreases very rapidly as a function of 
777,55 foi' values of 77755 = 6.5 x 10^^ (4.5 x 10^^) pmSugra and type-II can be formally 
distinguished by more than 3 (5) standard deviations. 

Also the errors A(7r7o) and A(Mi/2) do show dependence on 77755, especially at larger 
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FIG. 8: Error of mss, against 77155 for 5 free parameters. For these plots we used the LHC and 
ILC observables. The chosen values of parameters correspond to MSP-2 (left) and MSP-3 (right). 
The plots show the results for Seesaw type-II. 



values of mss- Again, the reason for this dependence is found in the strong correlation 
among those three parameters, as discussed above. The error A(tan /3) (and to some extend 
A{Ao)), on the other hand, shows less dependence on 77155. This is explained by the fact 
that the lightest Higgs mass, rriho, shows very little dependence on the seesaw scale. The 
slight dependence of A{Aq) on 77755 is due to the lightest stop mass. For simplicity in all 
following plots we show only the 1 a allowed regions. 

Fig. ( |8]) shows two more examples of A (77755) as a function of mss- Here results for 
the points MSP-2 and MSP-3 are shown for type-II seesaw. Only A (77755) as a function 
of 77155 is shown. We do not repeat the plots for the other parameters because they are 
qualitatively very similar to the case shown in fig. (^^. As the plots show results for MSP- 
2 and MSP-3 are similar to MSP-1. Values of 77755 below roughly 77155 ~ 10^^ GeV are 
inconsistent with pmSugra. This implies that for the ILC errors as estimated in j3] and jsl a 
combined ILC+LHC analysis should be able to distinguish pmSugra from type-II seesaw for 
nearly all values of 77755 relevant for neutrino data. We stress that this conclusion is correct 
only for those mSugra parameters for which (both left- and right-) sleptons and the lightest 
neutralino are kinematically accessible at the ILC. 

Up to now we have shown only results for seesaw type-II. Fig. shows a corresponding 
calculation for type-Ill and mSugra parameters as in MSP-1. Again, MSP-2 and MSP-3 show 
similar behaviour and we do not repeat the plots for these points. Again, the scale of 77755 is 
different from the case of type-II. Since SUSY masses show a stronger dependence on 77755 
in type-Ill than in type-II, larger values of 77755 can be distinguished from pmSugra in this 
case. In the examples shown in the figure all values of 77755 below roughly 77755 ~5 — 6xl0^^ 
GeV can be distinguished from pmSugra with more than 1 a c.l. Recall that in type-Ill one 
expects 77755 ^ 8 X 10^^ in order to explain neutrino data. Such "low" values of 77755 differ 
from pmSugra in the fits by many standard deviations. 

Errors on mo are similar to the values observed for type-II, while A(Mi/2) is larger in 
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FIG. 9: Error of mss^ 't^o and against for 5 parameters varied freely. For these plots 

we used the LHC and ILC observables. The chosen values for mss, itt-q, ^1/21 tan(3 and ^0 are 
the values according to MSP-1. The plots show the results for seesaw type III. 



type-Ill. Tlie correlation between mo, M1/2 and mss, discussed above for type-II, is also 
present in type-III and with an even larger correlation between M1/2 and mss in this case. 
The mSUGRA solution does not show up explicitly as a second separate minimum, but 
deforms the distributions, thus cutting the allowed ranges of mo and M1/2. 

Up to now we have always used a seesaw spectrum as input. One can also ask the 
opposite question: Can a pmSugra point mimic a seesaw spectrum? An example of such a 
calculation is shown in fig. fllUp . In this figure we show the allowed ranges for mss and mo 
for mSugra parameters as in SPSla' for type-II (left) and type-III (right). As one can see 
mss as low as mss ~ 10^^ GeV {mss ~ 7 x 10^^ GeV) are allowed at 1 a c.l. for type-II 
(type-III) fits. Also note that A(mo) is much larger than in a pmSugra fit, due again to the 
observed correlation among parameters. The results shown in fig. (fTOj) are consistent with 
the results discussed above, when a seesaw spectrum is used as input: mss compatible with 
Mq is reached at a very similar value of mss- 

Finally we note, that distinguishing type-II from type-III requires extremely high pre- 
cision, since they differ only at 2-loop order. The reason is that for 1-loop RGEs one can 
always cancel the shifts in the coefficient of the beta-functions by a rescaling of mss- We 
have checked this numerically. 

Closing this section we note that all results shown above have been obtained for the full 
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FIG. 10: The plots show random walks in which as starting point SPSla' was chosen. For the 
parameter fit we used mSUGRA plus seesaw type-II and III, respectively. The runs take into 
account LHC and ILC observables. 

2-loop calculation. We have repeated the exercise in several cases using 1-loop RGEs only. 
As a general result, due to the weaker running of 1-loop RGEs, differences between pmSugra 
and seesaw are slightly smaller, leading to slightly larger errors on the parameters. For the 
case of the ILC+LHC analysis, however, differences between both calculations are rather 
small, with errors in parameters typically increasing in the order of (10-30) % when going 
from a 2-loop to a 1-loop calculation. 

D. LHC only 

In this subsection we discuss the results for an analysis using only LHC measurements. 
At the LHC observables do not measure SUSY masses directly. Instead, observables measure 
either mass differences or, in case of the edge variables, combinations of mass squared dif- 
ferences. Also one expects that LHC measurements will be much less precise than what can 
be done in case of the ILC. As a result the distributions for an LHC-only analysis show 
more complicated features than for the case of LHC+ILC. Especially it should be noted that 
in some cases we do not have a sufficiently large number of well determined observables and 
fake minima can appear, which will lead to sometimes rather large error bars on parameters, 
as discussed below. 

Fig. (ITT!) shows error bars on m^^, tuq, M1/2 against mss for the point MSP-3 and seesaw 
type-II, again for all 5 parameters varied freely. Note the change in the scale for mss, the 
largest value shown is mss = 10^^ GeV. For larger values of mss type-H seesaw can no longer 
be distinguished in this fit from pmSugra with at least 1 a c.l. Note, however, that A{mss) 
decreases very rapidly for decreasing values of m^^ and for values of m^^ below mss ^ (few) 
10^^ GeV pmSugra and type-II are formally different by several standard deviations. 

The figure shows also that A(mo) and A(Mi/2) are much larger for the case of using only 
LHC observables than in the combined ILC-I-LHC analysis, as expected. Errors on mo and 
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varied. For these plots we used only the LHC observables. The chosen values for 77155, 77io, -M1/2, 
tanii and Aq are the values according to MSP-3. The plots show the results for seesaw type II. 



Ml/2 decrease in general with decreasing mss- The increase in A(mo) and A(Mi/2) around 
fT^ss ~ 7 X 10^^ GeV is due to the appearance of a fal<;e side-minimum. Such falie minima 
appear only for certain ranges of 7/755. Depending on which side of the real minimum they 
appear they can lead to an asymmetric increase of the errors as observed in this figure. 

Fig. ( IT2I) shows an example of a corresponding fit for type-Ill. Again, A (7/70) and A(Mi/2) 
and A(7?755) are shown as a function of 77155 for MSP-3. Note that other points show qualita- 
tively similar behaviour and that the range shown for mss is comparatively small. Values of 
mss larger than mss ~ 5 x 10^^ are 1 a c.l. consistent with pmSugra. Since 7^55 < 8 x 10^^ 
to explain neutrino data, LHC-only can probe interesting parts of the parameter space, but 
certainly will not be able to cover all possible values of 77755 - unless LHC errors on mass 
measurements can be improved compared to expectations by considerable factors. 

For decreasing mss errors again decrease in general. There are two exceptions from this 
general rule in this plot. First, errors increase around 77755 ~ 2 x 10^^ GeV. This is again 
due to the appearance of a iake side minimum, which slowly disappears again when going 
towards smaller values of 77755. The large increase in the error bars around 777.55 ~ 6 x 10^^ 
GeV is due to the fact that for smaller values of mss in this calculation x% is lighter than 
e/j, i.e. the edges variables are lost completely. With only a few observables in the fit, all 
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FIG. 12: Error of msst "^0; -^1/2 against mss for 5 parameters varied freely. For tliese plots we 
used only the LHC observables. The chosen values for mss, "n^o, ^1/21 tan/S and Aq are the values 
according to MSP-3. The plots show the results for seesaw type III. 



based on mass differences, tuq and M1/2 can hardly be fixed at all. The dramatic increase 
of the error bars of mo can be understood easily from Eq. ( IT5l) . As this equation shows the 
sfermion masses behave approximately like = tuq + aM^i^- When all edges are lost the 
remaining LHC observables can be fitted by varying 77255 and M1/2 only. 

Finally we have calculated the allowed parameter space in a seesaw fit when the true 
input point is pmSugra. Two examples are shown in Fig. (fT3ll . The mSugra parameters are 
for SPSla' and type-II (type-Ill) seesaw is shown to the left (right). The allowed regions are 
much larger than in the combined ILC+LHC analysis, compared to the discussion in the last 
subsections. For type-II (type-Ill) values of mss as low as 777,55 ^ 10^'' GeV (77755 ^ 6 x 10^^ 
GeV) are allowed at the 1 a level. This is similar - and consistent - with the results discussed 
above for the opposite fit. 

In summary mass measurements from the LHC only should be able to distinguish between 
pmSugra and type-II (type-Ill) seesaw for seesaw scales below roughly 77755 ^ 10^^ GeV 
{jnss ^ 6 X 10^'' GeV). This conclusion depends critically on the possibility to measure 
accurately several observables, as we are going to discuss next. 
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FIG. 13: The plots show the result of a random walk in which as a starting point SPSla' was 
chosen. For the parameter fit we used mSUGRA plus seesaw type II and III, respectively. The 
runs were made for LHC observables enabled only. 
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FIG. 14: Error of mss with respect to the error of the ILC observables, fiLC- 



E. Required accuracies on mass measurements and A (77155) 

All our results shown above crucially depend on the size of the expected error bars for 
the different observables. In this section we therefore discuss in some detail: (a) Which are 
the most important observables in our fits? And, (b) How accurately do we need to measure 
them to distinguish pmSugra from seesaw for a given, fixed value of mss- Again we will 
discuss the ILC+LHC case first. 

Fig. f lT^ shows A(m.55) for the points MSP-1 for two choices of mss- A(m55) is shown 
as a function of the error of the ILC mass measurements. According to |8| it is expected that 
the ILC can measure SUSY masses of cind / states kinematically accessible with errors of 
the order (0.5-2) per-mille. We define a common factor //lc and multiply all relative errors 
given in table 6 of [sl with this common factor. A (77155) is then shown as a function of this 
factor in fig. (fT^ . Note that in this calculation we keep all LHC errors unchanged at their 
"standard values" . 

As fig. (JH]) to the left shows, A(?7i55) increases with the assumed errors of the ILC 
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FIG. 15: In this plot the errors of niss with respect to mss are shown. The different hnes belong 
to different runs where different combinations of ILC observables were switched off. The different 
lines belong in each case to the same mss values but were a bit to be able to distinguish them. 
The real mss value is the value of the pink error bars which correspond to the run for all ILC 
observables enabled. 



measurements. However, for MSP-1 and mss = 10 GeV, LHC measurements alone are 
sufficient to distinguish type-II from pmSugra. Thus, error bars on mss hardly increase 
going from fjLc = 24 to fuc = 32. This means that ILC data dominate the fit until errors 
are about one order of magnitude larger than estimated in jsl, for larger ILC errors LHC 
measurements become more important for this choice of mss- 

The situation is quite different for mss = 2 x 10^^ GeV, see fig. f ll4p right. While 
fiLc = 6 still allows to distinguish between pmSugra and type-II, for fuc = 8, A(mss) 
becomes to large to differentiate between type-II and pmSugra. The required accuracy of 
measurements of SUSY masses at the ILC is therefore a strong function of mss itself. Errors 
of the order (1-2) percent are in general tolerable for seesaw scales below mss = 10^^ GeV, 
while per-mille level measurements are required in the interval [10^^, 10^^] GeV. We note 
that other SUSY points behave very similar and that for type-Ill correspondingly larger 
errors are tolerable. 

Fig. f lT^ treats all ILC observables equally. An interesting question to ask is, of course, 
which ILC observables are the most important ones for our analysis. Fig. (fT5|) provides the 
answer. Again for the point MSP-1 and for seesaw type-II we show A{mss) as a function 
of mss for different calculations taking into account different observables. We have kept all 
LHC observables "on" at their standard errors. "All ILC obs" is the standard fit, taking 
into account all kinematically accessible mass measurements with their original errors from 
[sl. We then switched off by hand completely the contributions from different observables. 
Switching off the measurement of the mass of ti hardly changes the result. On the other 
hand, it can be seen that measuring left-slepton masses is highly important. Error bars 
increase sizeably if this observable is not taken into account and while a set of measurements 
with all observables can distinguish pmSugra from type-II all the way up to mss = 10^^ GeV, 
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FIG. 16: Error of mss with respect to the error of {nig — mj^^) and (m;/)''*^^^. According to [7[| the 
error on {nig — m^^) is expected to be ~ 3.5% and the error of {muy^^^ is estimated as ~ 0.17%. 
Note that we also changed the error of {nig — m^_^ ) , accordingly. 



without the accurate measurement of mj^ all values of 7715-5 > 4 x 10^^ GeV are compatible 
with pmSugra at the 1 a level. The relative importance of II despite its larger error stems 
from the fact that mj^ has very little sensitivity to 77155, s^e fig. ([T]). 

We now turn to the discussion of LHC errors. In this case we do not use any input from 
the ILC. Fig. ( !T6|) shows A(7?755) as a function of the assumed error. Two observables are 
shown: to the left as a function of A-^ and to the right as a function of the edge variables. 
Note that for SPSla the LHC error for A-^^ is estimated to be ~ 3.5 %, while A(77i^;*^^'^) 
should be measured to an accuracy of 0.17 %. Note that, while Aim^p") can possibly be 
accurately measured in wide ranges of mSugra parameter space, the accuracy with which 
A-^ can be measured is far less certain. Both smaller and much larger errors on this quantity 
have been found in different study points, see 0]. 

The figure shows that for mss = 2 X 10^3 a 7 % error on A-j^ (which corresponds to a 
factor of 2 in the plot) is sufficient to distinguish between pmSugra and type-II, while an 
error of 9 % on this quantity is not sufficient. Again, the maximum value of this error which 
still allows to distinguish between type-II and pmSugra is a strong function of the (unknown) 
77755 itself. However, we have found that always A-^^ is a critical input observable for our 
analysis. ^ The importance of A-^^ can be understood from Fig. ([2]) and (|3]): coloured 
sparticle masses depend much more strongly on m55 than masses of, for example, sleptons. 
Thus, despite the larger relative error on A-^;^ compared to the edge variables, it is nearly as 
important as demonstrated in fig. ( |T6l) to the right. 

Fig. (flTll shows the results of different runs, where we have switched off artificially different 
combinations of observables. As noted above, A-^^ and the edges are the most important 
observables for fixing A (77155). However, the Higgs mass measurement is not negligible, 
despite the fact that A (77755) does not increase much in the figure, when rriho is switched off. 



We also consider A-^^, which, however, is less iniportant due to its larger error. 
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FIG. 17: In this plot the errors of 77155 with respect to mss are shown. The different hnes belong 
to different runs at which we switched off different combinations of LHC observables. The different 
lines belong in each case to the same mss values but we separated them in the plot a bit to be 
able to distinguish them. The real mss value is the value of the pink error bars that correspond 
to the run in which the Higgs mass measurement is disabled. 
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FIG. 18: In this plot the errors of mss with respect to mss are shown. The different lines belong 
to two different runs at which we set the Yukawa couplings to very small values (~ 10~^) in the 
first run and in the second one we fitted the Yukawa couplings to neutrino data and calculated the 
error with the fitted values. 



This importance lies in the fact that without m/jO the largest value of 77155 not compatible 
with Mg is 2 X 10^^ GeV, compared to 77755 = 10^^ GeV for rriho switched on. 

We do not repeat the discussion for type-Ill seesaw. Results are very similar qualitatively, 
but again larger values of 77755 can be tested for the same errors on the observables. 

Finally, we turn to the question of Yukawa couplings. Fig. ( fTSj) shows again A(7?755) as a 
function of mss for two different calculations: (i) a calculation with triplet Yukawa couplings 
negligibly small (all (It)*^ ~ O(10~^)) and (ii) a calculation in which Yr has been fitted to 
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give the atmospheric and solar neutrino mass squareds at their best fit values with neutrino 
angles taking tri-bimaximal values. As can be seen, differences between both calculations 
become negligible below roughly mss = 10^^ GeV, as expected. For larger values of mss 
correctly fitting the Yukawas leads to slightly smaller errors on mss- This can be understood, 
since for finite Yukawas SUSY masses change slightly stronger than for infinitesimal values 
of Yt, making the fit easier. Note, however, that we have not scanned over all allowed range 
of Yt in this calculation. In a complete 14 parameter fit errors might be larger. Note 
also, we can not find any good neutrino solution for rriss larger than rriss = 6 x 10^^, since in 
this calculation we have chosen for the coupling A2 = 0.5. In conclusion, a full fit including 
Yukawa couplings will be necessary only if signs of mss ^ 10^"^ GeV have been found in 
SUSY mass data. 

IV. CONCLUSIONS AND DISCUSSION 

We have studied the possibility to obtain indirect information on the seesaw scale from 
SUSY mass measurements at future colliders. Since in the type-I seesaw only SM singlets 
are added to the MSSM particle content, none of the measurements which we considered are 
expected to show sizeable departures from mSugra expectations. We therefore concentrated 
on the type-II and type-Ill realizations of the seesaw. 

Assuming mSugra boundary conditions and taking error estimates as forecasted by study 
groups we find that a combination of LHC and ILC measurements should be able to distin- 
guish pure mSugra from mSugra plus either type-II or type-Ill seesaw for nearly any relevant 
values of the seesaw scale, if (a) at least Xi? (^bJ P'R and e^/ fiL are kinemetically accessible 
at the ILC and (b) the LHC can measure mg — m^^ and the edge observables accurately. We 
always assume that the lightest Higgs has been found. The degree of confidence with which 
pmSugra can be distinguished from mSugra plus seesaw depends sensitively on the actual 
value of the seesaw scale. At the "critical" value of mss ~ 10^^, beyond which neutrino data 
can no longer be explained with Yukawas smaller than 1, the difference between type-II + 
mSugra and pmSugra could be as low as only 1 a c.l. However, the difference between pm- 
Sugra and mSugra -|- seesaw rises very sharply with decreasing 777,55 and formally more than 
5 cr c.l. could be reached already at (5 — 6) x 10^"^ for type-II. Differences between pmSugra 
and mSugra plus type-Ill are always found to be larger than for mSugra plus type-II for the 
same value of m^^. 

As expected, the future is not as bright, if we take into account only LHC data. Nev- 
ertheless, using only LHC data one can distinguish pure mSugra and mSugra plus seesaw 
in some favourable parts of parameter space. Especially, we point out that the lower the 
real value of the seesaw scale mss is, the easier it becomes to distinguish pure mSugra from 
mSugra plus seesaw. We have discussed the most important measurements for the LHC and 
the ILC and the relative errors with which these observables need to be measured for this 
analysis to be possible. In our analysis we used exclusively mSugra SUSY breaking bound- 
ary conditions, but other, more involved SUSY breaking schemes with more free parameters 
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could, in principle, be "tested" in a similar way. 

Of course, the analysis presented in this paper is far from being complete. If SUSY is 
found at the LHC, one would need to redo our calculations with real data. However, the 
experimentalist would not, as we have always assumed in our fits, know the real values of 
the parameters. We have tried for a few points, whether the correct input parameters can be 
retrieved for arbitrary starting points in our MC random walk procedure and - given enough 
CPU-time - are able to find the correct minimum. However, our "observables" are theoreti- 
cally calculated observables and thus perfect in contrast to real data which are expected to 
scatter around the true values and might show tension between different observables. Thus, 
finding the correct minimum in real data might be more difficult. Moreover, the underlying 
model will not be known a priori and, thus, xl/m different models need to be calculated 
and compared. 

Nevertheless, even taking into account the limitations of our study, we think it is highly 
motivating that type-II and type-Ill seesaw leave sizeable traces in SUSY spectra, which 
should show up, if sufficiently accurate mass measurements are possible and become avail- 
able. 
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